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Abstract 

Generic 50 (10) GUT models suffer from the problem that Planck scale induced non-renormalizable 
proton decay operators require extreme suppression of their couplings to be compatible with present exper- 
imental upper limits. One way to resolve this problem is to supplement 50 (10) by simple gauged discrete 
symmetries which can also simultaneously suppress the renormalizable R-parity violating ones when they 
occur and make the theory "more natural". Here we discuss the phenomenological viability of such models. 
We first show that for both classes of models, e.g the ones that use 16h or 126// to break B-L symmetry, 
the minimal Higgs content which is sufficient for proton decay suppression is inadequate for explaining 
fermion masses despite the presence of all apparently needed couplings. We then present an extended 16h 
model, with three 10 and three 45-Higgs, where is free of this problem. We propose this as a realistic and 
"natural" model for fermion unification and discuss the phenomenology of this model e.g. its predictions 
for neutrino mixings and lepton flavor violation. 
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I. INTRODUCTION 



The neutrino observations of the past decade have put the spotlight on gauged B-L symmetry 
as well as unification groups such as SO(IO) and SU{2)l x SU{2)ji x SU{A)4 containing B-L 
as prime candidates for theory of matter, forces and flavor. While both these groups incorporate 
the seesaw mechanism for neutrino masses, SO(IO) has the additional attractive feature that gauge 
couplings unify at high scale. It is however highly nontrivial to obtain a "truly natural" SO(IO) 
model due to such issues as doublet triplet splitting, rapid proton decay etc. In this paper we 
discuss how one aspect of this naturalness can be addressed i.e. how one can naturally suppress 
proton decay in SO(IO) models while preserving our understanding neutrino masses. 

We first note that SO(IO) models for neutrinos discussed in recent literature can by and large 
be divided into two classes: 

(i) One class which uses only renormalizable couplings involving the Higgs fields 10, 120 and 
126 for fermion masses and the last multiplet for breaking B-L symmetry and multiplets such as 
45 and/or 210 for gauge symmetry breaking yj]. This theory could be considered as an ultraviolet 
complete theory by itself. 

(ii) The second class uses 10 plus 16 © 16 for fermion masses with the 16's breaking the B-L 
symmetry. Here one generally uses 45-1-54 Higgs fields for SO(IO) breaking. An important feature 
of this class is that it has to rely on nonrenormalizable couplings to understand fermion masses 
and therefore has to be viewed necessarily as an effective theory at the GUT scale[|20. 

The first class of models leads to automatic R-parity conservation when SO(IO) breaks down 
to MSSM so that there is a natural candidate for dark matter whereas the second class of models 
suffers from R-parity breaking and hence has no stable dark matter in the absence of additional 
symmetries. So in principle one could argue that this class of models are not "pure" SO(IO) 
models. 

Both models have an additional naturalness problem arising from the fact that they allow R- 
parity conserving nonrenormalizable couplings of the form X16^/ Mpi which lead to rapid proton 
decay. Such interactions could be induced by nonperturbative Planck scale effects and it is there- 
fore not safe to ignore them. Present proton life time limits constrain A to be < 10^. Such a 
small value of A would suggest that there is probably a symmetry responsible for its smallness. 
This question is particularly urgent for the class of SO(IO) models with 16 Higgs since they rely 
on other such dimension four higher dimensional operators with coefficients of order one to un- 
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derstand fermion masses. This problem is generic to all non-GUT susy theories such as MSSM 
or left-right models as well as SU{2)l x SU{2)r x S'f/(4)4 models and not just GUT theories. 
One way to understand the suppression of such operators despite the presence of non-perturbative 
gravitational effects, is to have an additional gauge symmetry beyond SO(IO) which can forbid 
these unwanted terms. The simplest possibility is to have a discrete gauge symmetryyO. There are 
of course other possibilities [|4|]. 

The discrete gauge symmetry supplemented SO(10) models that suppress proton decay were 
studied for a large class of models in a recent paper[|5D. In particular two minimal SO(IO) models- 
one with 16-Higgs breaking the B-L symmetry and another with 126 breaking B-L were shown to 
be free of both proton decay problem as well as R-parity problem if SO(IO) was supplemented by 
a gauged Zq symmetry. They looked promising for phenomenology since all necessary terms in 
the superpotential for phenomenology were allowed by the symmetry. It is the goal of this paper 
to study the viability of these models. 

The results of this paper are the following: (i) the minimal versions of both IG^-based as well 
as 126-based models discussed in Ref.ljsl] are not realistic since they fail to give desired MSSM 
doublets that would be required to give rise to realistic fermion masses and mixings; (ii) if the 
16-based models are extended to have three 10-Higgs fields and three 45 multiplets, one can have 
the desired doublet-triplet splitting and fermion masses that can match observations. This model 
differs from other 16-based models in that proton decay here arises only from the gauge boson 
exchanges unlike other models where Planck scale induced effects as well as Higgsino exchange 
ones play a rolej^j]; (iii) we study the phenomenological implications of this model and isolate 
some of its tests e.g. in the domain of lepton flavor violation. 

This paper is organized as follows: in sec. 2, we review the salient features of the two classes 
of models; in sec. 3, we discuss doublet-triplet splitting problem of the minimal models; in sec.4 
we discuss the three Higgs extension of the 16-based model that fits fermion masses and mixings; 
in sec. 5, we discuss how large neutrino mixings and observed neutrino masses arise in this model. 
We summarize our results in sec. 6. 

II. THE S^O(IO) X MODEL FOR 16-HIGGS B-L BREAKING 

The main features of generic SO(IO) models with 16-Higgs fields breaking B-L symmetry are 
the following: (i) the quarks and leptons are assigned to three 16-dimensional spinors (denoted by 
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ijm, m=l,2,3); (ii) the GUT symmetry is broken down to SU{2) l x SU{2)r x f/(l) b-l x SU{3)c 
by a 45 © 54 set of Higgs fields; (iii) SU {2)r x f/(l) b-l symmetry is broken by the 16-Higgs pair 
denoted by ipn © tpH set to the standard model symmetry which is then broken by SU (2) doublets 
that are linear combinations of the doublets in SO(IO) 10 and 16-Higgs fields. The standard model 
symmetry along with supersymmetry emerge just below the GUT sclae of 2 x 10^^ GeV. 

This model has several naturalness problems: it not only allows the dangerous R-parity con- 
serving ^^g^ terms but also terms such as '''"'^X/p^^ terms which on B-L symmetry breaking lead 
to all three types of R-parity violating operators present in general MSSM i.e. LLe'^, QLd^ and 
u^df^d^ type. Thus this model for natural values of couplings will lead to extremely rapid proton 
decay which is unacceptable. The question addressed in Ref.dsj] is to search for gauged discrete 
symmetries that will keep the model phenomenologically viable while keeping them "proton decay 
safe" and it was shown that the minimal anomaly free discrete gauge symmetry is Z^. Similar con- 
siderations for 126 type models also led to the symmetry and in both cases an extra 10-Higgs 
field denoted by H' in addition to those considered already. 

To see the discrete symmetry charges for various fields that forbid both R-parity violating terms 
as well as R-parity conserving baryon number violating terms, while at the same time keeping the 
required terms responsible for good phenomenology, we divide the superpotential terms into two 
classes: type I terms that must be kept for phenomenology and type II terms that must be forbidden 
to suppress proton decay and R-parity violating terms. They are given below: 
Terms of type I: They include il^m ipm H, {ipm ipnY /^^^ i^n "^h^ S"^'^ ^nd S A^, where H, A, 
S are 10-, 45-, 54-plets, respectively. Taking the discrete gauge symmetry to be Zn, we can write 
down the constraints on the Zn charges that are required by the type I terms: 



2q^^ +qH = 





mod , 


2g^„ + 2q^^ = 





mod A^ , 


(la) 







mod , 


Qh + qw = 





mod A^ , 


(lb) 


2qA = 





mod , 


2qs = 





mod A^ , 


(Ic) 


3qs = 





mod , 


2qA + qs = 





mod A^ . 


(Id) 



Here, we denote the Z^ charge for a field Fhy qp- 

Type II terms: These are the terms that must be forbidden from appearing in the superpotential 

and are ip^ tp^ H', ip^, 'ipn and ipl^ ipn, 'ipm i'H H, ipn H' and ip^ip^A. We forbid the 
ipm 1pm H' in order to avoid large Higgsino mediated contribution to proton decay since this is 
the very problem we are trying to solve. The necessary constraints on the Zj^ charges have to be 
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chosen such that they satisfy the inequalities 



2g^™ + gH'^0 modiV, 4gv.„7^0 mod iV , (2a) 

%m + 'i^„l^ mod N , 3 g^-™ + ^ mod N , (2b) 

+ (lH',H + q^„^ mod , + q^^ + ^ mod . (2c) 

The last set of constraints come from the requirement that the discrete symmetry must be a 
gauge symmetry i.e. it must be anomaly free. The anomaly freedom constraints are: 

16 {Ng q^^ + q^^ + (t^^) + 10 {qu + qw) + 45 gA + 54 = mod N' (3a) 
2A^g + 2 g^,^ + 2 g^^ + gn + qw + 8gA + 12g5 = mod (3b) 

{N, odd N 
(3c) 
N/2, even A^ 

It was shown in [SI] that the smallest symmetry allowing us to fulfill all criteria is for number 
of generation Ng = "i. A possible charge assignments is g^„ = 1, g^^ = —2, q^^ = +2, 
qH = —2, qw = +2, g45,54 = (cf. tables |I](a) and (b)). This charge assignment allows for seesaw 
couplings and the possibility of fermion masses from couplings of type ^l^m H. The allowed 
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operator ipm "ipm i^u contributes to both the fermion masses as well as to the seesaw. The model 
also eliminates the dangerous proton decay operator QQQ Lor operator of type {'ijjm)'^/Mp. 

While the allowed set of operators provide a necessary condition for the model being phe- 
nomenologically viable, the final step where we judge whether it is acceptable first requires that 
we do doublet triplet splitting and see if the sub-GUT scale structure of the model can generate 
acceptable pattern of fermion masses or not. We address this question in the next sub-section. 

A. Phenomenological Viability of the 16-Higgs model 

To analyze the phenomenological implications of the model, let us start by writing down the 
superpotential allowed by the discrete symmetry and SO(IO) invariance: 

W = Wy + Wh 

Mpi Mpi 
Wh = MHH' + SHH' + AHH' + ijH^HH + i)Hi>HH' + M^i)^ (4) 



where a denotes the various irreducible representations in the product of two 16's. 



TABLE I: 



(a)MSSM part (b)16-ffiggs model. (c)126-ffiggs model. 



Field 


quantum numbers 




Field 


quantum numbers 




Field 


quantum numbers 




16i 






16_2 




A 


1262 


H 


10_2 




i'H 


T62 




A 


126 2 


H' 


IO2 




A 


45o 




E 


2IO0 






S 


54o 







First we want to find whether we can solve doublet -triplet splitting problem within this model 
field content. The vev s of the 54 and 45 can be assumed to have the following forms: 



<A>-- 



Diag{a, a, a, b, h) 



3 3 

< 5 >= ( \ ) Diag{s, s, s, --s, --s) 



It is useful to express all fields in terms of the SU(5) multiplets. 



(5) 



16= 1 + 5 + 10 
16 = 1 + 5 + To 
10 = 5 + 5 



(6) 



so the mass matrix in terms of SU(5) components of the fields will look like 



(5ij,5ij',5:^^) 



M+A+5 c 
M-A+S 



•^11 







5h' 



(J) 



where c is the vev of the 16 and 16 of the SO(IO), c = (V'i?) = (V'i?)- ^iid to get expression for 
the mass matrix of the doublets(triplets) one has to substitute instead of A(S) for doublets 

and a{s) for tripets respectively. One can see that this matrix can have a zero eigenvalue only if its 
determinant vanishes i.e. 



Det = (M - A + S) (c^ - M^(M + A + ,5)) = 0; 



(8) 
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This equation has two solutions: taking the first one i.e. M — ^4 + 5" = 0, we find for the doublet 
mass matrix Mud in the basis {H,H' , t/jn) to be; 



2b 



^ z c ^ 

c M 



2{M--s) 



(9) 



To find the MSSM doublets in terms of the GUT submultiplets, we diagonalize M^d and find 
its zero mode eigen-vector. The usual MSSM Higgs fields /i„ will be linear combinations of 
{H, H' , ipff) that correspond to the zero mode eigen-vector of the above matrix. From the follow- 
ing equations, we find: 



\ 



z c 

^0 c M 

^ 

z c 
c M 



^ D ^ 

D21 

U21 



(10) 



It is easy to see that 



hd — Hd 



(11) 



Since the 10-Higgs denoted by H' does not couple to 16-fermions, the MSSM up-Higgs doublet in 
this model does not couple to matter and therefore all the up quarks remain massless. The solution 
(i) to the determinant equation is therefore not acceptable. 



Turning now to the second solution i.e. (ii) M + A + S 
mass matrix: 

\ 

Mud^ 



M, 



-, we get the structure of the Mud 



X 
c M 



(12) 



In this case the zero mode corresponding to the MSSM doublet hd, can be represented by the 
column vector: 

/O \ 



D 



M 



/ M 



u 



m2 







(13) 



The MSSM doublet hd in this case does not couple to quarks and charged leptons and leave those 
fields massless. Again this is not acceptable. Taking these two cases together we conclude that in 
the minimal gauge discrete symmetric Zq model, the doublet triplet splitting and nontrivial fermion 
masses cannot happen simultaneously and the model is therefore not phenomenologically viable. 
We wish to emphasize again that in drawing this conclusion, we have also considered higher 
dimensional operators that could contribute fermion masses. It turns out that in our, case operators 
such as {ipmYiipHY, ii^mYH'A which can lead to the bottom quark mass are not allowed due to 
the Zq charge assignments. 



B. Phenomenological viability of the 5*0(10) x Zq 126 model 

This class of models typically have the Higgs multiplets of type 10 126,210 (and 120) fields 
to explain fermion masses including neutrino masses and mixings. These models do not have R- 
parity breaking terms even after GUT symmetry breaking. The Zq charge assignments that makes 
the model proton-decay- safe while keeping necessary terms for possible fermion masse are given 
in Table 111(c). 

The Higgs superpotential in this case looks like 



One might think that this can lead to a realistic model for fermion masses. However as in the 
previous sub-section, we must analyze the doublet-triplet splitting in order to study the fermion 
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masses, we will see that in this case too there is a conflict between the doublet-triplet splitting and 
fermion masses. 

To see this, we write down the mass matrix for MSSM doublets contained in various GUT 
Higgs multiplets in the theory: (see ||7|] for the exact Clebsch -Gordon coefficients) 



i^Huy H^y 

) X Mud X 



H'd 
A, 
A, 



ud 



1 





M„ 







^HAE^. \ 




Mh 





^(*-^) 




















I 












M. + ^($, + ^) 


■^AAE ,. 
10 


V 















15) 



Where $i, $2, '^'s are the vevs of S in different directions and f , f vevs of A, A fields respectively. 
For simplicity of analysis, we rewrite this matrix in the symbolic form as follows: 

/O m ^ c 

b 

m b 

a M 

d M2 A 

\ ci A B J 
We want to have one massless state so we require that 



(16) 



DetM, 



ud 



0. 



(17) 



This requires that one of the following two conditions be satisfied a = or a = 

Abcrd-AHrn-bcc^Nh+BbmM2 _ ^ow investigate the first one (i)a = ^ which leads to the 

following zero mass eigenstate: 

/o \ 



u 



M 



b^l+ 
1 



'1+ 



Al2 








(18) 



D 



(,2 



(19) 



V 



In this case we see that the up-MSSM Higgs doublet does not couple to matter fermions. 

Turning now to the case (ii) where we have a = Ahc^d-AHm-bcM+BbmM^ ^ f^j. same 

eigenstates: 

/*\ 



U = (20) 

* 

V*/ 

* 

D = (21) 

V*/ 

where we have written only the zero entries in the columns. The *'s represent non-zero entries 
whose detailed form is irrelevant for our discussion. It is clear that in both cases the doublet -triplet 
splitting is incompatible with giving masses to the fermions; in the case (i) to up quarks and in 
case (ii) to the down quarks. 

So neither the minimal 16 nor 126 models when made proton decay safe can lead to viable 
fermion masses along with doublet triplet splitting. We therefore have to extend the Higgs sector 
to get a realistic model. In the next section, we give one such example and analyze its flavor 
phenomenology. 



m. EXTENDED 16-HIGGS MODEL 



We now extend the 16 -Higgs model by adding extra Higgs multiplets in such a way that 
anomaly freedom as well as proton decay constraints are satisfied and yet the model can lead to 
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viable phenomenology. The simplest possibility appears to be to extend our 16 model by adding 
one additional 10 - and two 45-^42, A3 fields under SO(IO), where has zero charge under 
Zq, and ^2,3 charges are -2 and 2 respectively (see table below). Zq charges are easily seen not to 



TABLE H: 



Field 


quantum numbers 


A2 


45_2 


A3 


452 


H3 


lOo 



ruin our anomaly cancellation conditions. We also have redefined our H', H fields as Hi^2 and A 
as Ai for the simplicity of notation . Now the superpotential is given by: 

W = MH^H2 + + A^H^H2 + ^h^hH^ + ^h^hH2 + M^i^ui^H + M3HI + 

+H2H3A2 + H.H^A^ + ""'^f" + (22) 

Our model allows an operator of the form '^^^'^ where substituting the vev of the field A3 we get 
a proton decay operator with effective A ~ <^ 1 but not suppressed enough to be acceptable. 
However this problem disappears if the vev {A3) — 0. We will see below that there is an allowed 
vacuum, where indeed this is possible. 

To study the doublet triplet splitting in this model, note that the mass matrix for the 5 and 5 of 
SU(5) is given by 

/ c\ 

M-Ai + S ^2 

-A2 M3 

c 5 M^y 



V 



Where 5 in the (43) element of the matrix comes from 



(23) 



ud 



H:iA2tpHipL 
Mpi 



coupling. As before, we want 



the determinant of this matrix to vanish. This leads to the following constraints 
Case (i) 



M-Ai + S = 



(24) 



Case(ii): 



M + Ai + S 



c^M3 



cAoS 



M3M, 



(25) 
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Here we consider only the simpler of the two cases above i.e. case (i) to illustrate that our proposal 
leads to a realistic model. In the first case M — Ai + S = 0,D = (1, 0, 0, 0) (implying that the 
hd has non-zero component in the multiplet Hi) in the same way as was in the minimal model, 
but now due to the presence of A2 field all the f/i is nonvanishing, so that the "up" quarks will get 
masses from the ^.^-f^i operator. In the next sections we will discuss the detailed fit to fermion 
masses for this extended 16 model. 

As we can see extended 16 model can solve doublet-triplet splitting problem as well as provide 
masses for all fermions, but now we have to check whether higgsino mediated proton decay oper- 
ators are allowed. Even though quarks and leptons couple only to the Hi field and there is no mass 
term oc HiHi, mixing between Hi, H2, H3, ipH fields can lead to the nonvanishing diagrams with 
higgsino exchange. The contribution of these diagrams will vanish if only the (HiHi) element of 
the inverse mass matrix (l23l) for the heavy triplets vanishes, thus the triplet part of the A2 should 
be zero. We will see in the next section that the requirement of the < ^3 >= combined with F 
flatness condition will lead to this condition. 

IV. SUPERSYMMETRY DOWN TO THE WEAK SCALE 

First we want to find out whether there is a minimum of the potential that can correspond to 
the solution we are interested in i.e. having supersymmetry survive down to the weak scale. The 
Higgs part of the superpotential is: 

Wh = M^ipH'ipH + miAl + 7712^2^3 + m,S^ + XiAjS + XiS^ + \3A2A3S + Xii^ni^nAi + A5A1A2426) 

-i-nonrenormalizable terms 

The vev of the 45, 54 and 16 fields will in general have the following form: 

(A) = (-1 ^) {ai,ai,ai,bi,bi); 

3 3 

{^Ph) = {i^H) = c (27) 
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So we can rewrite the superpotential in terms of the vev of these fields, using the further identities: 

miAl = -2{3al + 2bl)mi 

7712^3^2 = -2(3a2a3 + 26263)7712 
X^AlS = {-Gals + 66^s)Ai 

A3A2A35' = (-602035 + 662&3)A3S 
X2S' = -yA2.^ 

m^S^ = Ibrriss'^ (28) 
The condition of the vanishing F terms leads to the following constraints, 

— 12777iai — 12AiaiS + 3X40^ — 
-8777i6i + I2A161S + 2A4C^ = 
—6777203 — 6A3a3S = 

—47772^3 + 6A363S = 

—6777202 — 6A3a2S = 



9a 1 
dW 

'dh 
dW 

da2 
dW 

dW 
das 
dW 



db: 



-4777262 + 6A362S = 



'3 



^ = 30777,5 + Ai(66? - 6a?) + A3(-6a2a3 + 662^3) - y^2S^ = 
dW 

-— = 2cM^ + 2cA4(3ai + 26i) = (29) 
dc 

we are interested in whether there exist a solution with 03 = 63 = and 62 7^ these constraints 
lead to the following restrictions on the vevs 

27772 



02 = 0, s — 



3A3 



(3a, + 26OA4 = -M^, . = I'^f' (30) 

Ai(2ai + 36i) 

required to suppress higgsino exchange diagrams. Now we will present the other massless com- 
ponents of the higgs fields that provide the breaking of the SO{10) SU{3) x SU{2) x U{1). 
We will identify them by their charges under SU(3) x SU(2) x U(l) 
1) (3,1,2/3) fields (^1,^2,^3,^77); 
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/-4mi - 4Ais 

-2m2 - 2A3S + 2miA5 

— 2m2 — 2A3S — 2miA5 



-2A4C \ 



A4(26i-ai) + M^y 



(31) 



\^ -2A4C 

2) (3,2,-5/6) fields (^1,^2,^3,^) ; 

/ -4mi + Ais -ih2\-> 2z(ai + 6i)Ai\ 

-2m2 - A3I + i(ai - 6i)A5 

i62A5 -2m2 - A3I - i{ai - 6i)A5 zA362 

\—2i{ai + 6i)Ai —i\zb2 4:ms — 3A2S J 

3) (3,2,1/6) fields {A,, A^, A3, S, i^n); 



(32) 



/ — 4r?ii+Ais 


262^5 -2m2+A3|-i(ai+bi)A5 
-2i(ai~6i)Ai 
-2A4C 



4) (1,1,1) fields iAi,A2,As,ijHy, 



—162-^5 2j(ai— 6i)Ai — 2A4C \ 

-2m2+A3f+i(ai+fei)A5 * 

-i62A3 

ifeAa 4ms— 3A2S 

M^+A4ai / 



(33) 



/-4mi + 6A1S -2^562 

-2m2 + 3A3S + 2ibiX5 

2iX^b2 —2m2 + 3A3S + 2ibiX^ 



2A4C 





(34) 



\ 2A4C A4(3ai - 2bi) + M^J 

from the equations (I29H30I) one can see that each of these matrices will have one massless eigen- 
state. So we have total 32 massless goldstone bosons. One more goldstone boson needed to break 
SO(IO) down to SU(2)X SU(3)X U(l) comes from the phase of the , fields 



V. FERMION MASSES 



The following couplings allowed by Zq x 5*0(10) symmetries will lead to fermion masses after 
symmetry breaking. 



+ r 



PI 



5!M, 



PI 



3!M, 



PI 



3!m; 



PI 
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Where h^^, /^°, /^^^, g^^^ are symmetric 3x3 and fc^^^, (7^^° antisymmetric matrices. The upper 
index of /^^^..., 10, 120, 126 shows the SO(IO) structure of the fermion couplings. 



D 



M. 



Vu 



a — b 



M. 



+ 9 



120 



PI 



h{a - h) 



+ 9 



^^^ah + a^ \ , ^ (f^o^ _ ji26 



Ml, 



M. 



PI 



PI 



Vdl 



^10 _ ^^120*^ „120 



Mpi 
a + b 



120 b{3a + b) + 

~ 9 T7o, 9 



Ml, 



Ml, 




PI 



PI 



M, 



M 



3a — 6 



b{a + b) 
y » ^2 " 



-a6 + 

~Mir 



M. 



120 ^(3q - fe) i263(-a& + a^ 



Ml, 



Ml, 



16/126 V^_(^D)T(^g^l26 



pi 



M, 



PI 



(36) 



Where the vev of the fields and tpj^ are related to the vev of the MSSM doublets and 
in the following way 



(Hiu) = Oi{K) 



(37) 



Our claim is that these Yukawa coupling structure is rich enough to fit all the fermion masses. 
We give below an example of a scenario where correct fermion masses can arise. 

We take the case where all antisymmetric couplings vanish and that down quarks and leptons 
are brought to the diagonal basis at the same time. Thus g^'^^ and h^^ are diagonal, and all the 
mixing in the quark and lepton sector arise from the couplings /^^ and f^'^^. We now show that 
even under such limiting assumptions we can fit all the fermion masses. We know the quark 
masses at the GUT scale thus we can find corresponding h^'^, g^"^^, but on the other hand in the 
case when all the quark mass matrices are symmetric the mass matrix for the up quarks is equal to 



Mu = {CKMf.M^'''^.CKM; 



(38) 



this leads to the constraint on the linear combination of and /i^^.On the other hand we can 
find /i26 from the known neutrino masses and mixing angles so this fixes and /^^s. 
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So here is the fit for a,b,c,a,^ that leads to the good quark and lepton mixing. We set tan{f3) = 
55, the masses of the quarks and leptons and the vev's Vu, Vd at the GUT could be found in Isj]. 

a = 0.8 
e = -0.47 

7 = 1 

a = 0.028Mpl 
b = -O.OUMpl 
c = -0.024Mp/ 

(39) 



ho 

9l26 



f 



10 



126 



Dia^(0.000574402, 0.0208157, 0.792778) 
L)mc/(0. 114692, -4.43859, 0.190027) 

^-0.00858937 - 0.000536445? -0.0110172 + 0.00138491? 0.094032 - 0.039795i ^ 
-0.0110172 + 0.00138491 -0.286505 + 0.000708062i -0.508733 - 0.00956549? 
0.094032 - 0.039795Z -0.508733- 0.00956549? -0.134644 - 0.0176603i 



V 



-5.509 X 10-*^ + 5.992 x 10" 
0.00003070 - 0.00001463i 
-0.001095 + 0.0004482Z 



A 



0.00003070 - 0.00001463? -0.001095 + 0.0004482i 
-0.00014401 + 9.1768 x 10"^ 0.004929 - 0.0000360t4C|) 
0.004929 - 0.00003607? -0.1718 - 0.008830? 



these Yukawa couplings lead to the good mass matrices for the up, down quarks, charged leptons 
and neutrinos. Note that there are no dimension five operators that contribute to down quark mass 
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matrix due to the discrete symmetry of the model. 



M. 



CKM 



Dmc/(1.46323, 32.2949, 1638.17)MeV 
Dm^(0.35668, 75., 1636. )MeV 
Diag{0.757795, 208.466, 87029.8)MeV 
^ 14.74 - 3.341i -67.49 + 8.615i 586.2 - 247.9i 
-67.49 + 8.615 312.8 + 4.202z -3159. - 57.83z 
586.2 - 247.9i -3159. - 57.83i 86910. 
135016.MeV, Vd = 2065.81MeV 

Diag{9.577 * 10"^ 0.00007594, 0.002697)eV^ 
/ 



V' 



MeV 



-0.7995 0.59481 -0.083258 

-0.35036 - 0.14031i -0.51466 - 0.090935i -0.31226 + 0.69779z 
-0.46449 - 0.050427i -0.59179 - 0.15112i 0.2327 - 0.5954i 



0.973841 
-0.227079 - 0.000134603i 
0.00675874 - 0.00284968i 



A 



0.227198 0.00169092 - 0.00292876? 

0.97298 - 0.000031403? 0.0369876 (4lD 

-0.0364044 - 0.000664834i 0.99912 



The FIGlH shows the distribution for the values of sin^Or^ in our model for 2a values of neutrino 
masses and mixings. As is clear from this figure the model has a slight preference towards the 
region of small 6^13. The distribution for the other parameters of the neutrino mass matrices 
appear to be spread uniformly over the allowed regions. 

VI. CONSTRAINTS FROM THE LEPTON FLAVOR VIOLATION 



/ 



In this section, we discuss the predictions of this model for lepton flavor violation. As is well 



known 







Uj], even if the slepton mass matrices are diagonal at the GUT scale the ROE running 



down to the scale of the righthanded neutrino will lead to the mixing in the slepton sector, which 
via one loop diagrams leads to lepton violation. We will assume mSUGRA boundary condition 
for scalar partner masses and use the renormalization group equations to run them down to the 
seesaw scale when the right handed neutrinos decouple. 

We will work in the basis with diagonal righthanded majorana neutrino matrix, then the slepton 
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0.005 



0.01 



0.015 



0.02 



0.025 



sin" 01 



FIG. 1: Distribution plot for the values of sin^Oij, 



mixing will be approximately equal to 



LL 



(42) 



where Yy are the Yukawa couplings of the Dirac neutrino. These Yukawa couplings appear to be 
of roughly 

^ 10"^ 5 ■ 10"^ 5 ■ 10-3^ 
10-5 5 ■ 10-3 2 ■ 10-2 
10-5 3 ■ 10-3 0.3 



(43) 



Here Yy is a linear combination of the Yukawa couplings /iio, /lo and /i26 of the previous section. 
The slepton mixing leads to the lepton flavor violating processes U Ij'j with the amplitude equal 
to 



iM = emue^lj (iq^'axf.iALPL + ArPr)) k 



(44) 



Where the q is the momentum of the photon and Pl,b. = |(1 =F I5), the exact expression for the 
Al r can be found in [|ll|] . The branching ratio for this processes will be equal to 



487r'^a 



i\Al\ 



(45) 
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BR{fi 67) 


1.2 • 10"" 


BR{t n-f) 


6.8 • 10-8 



The present bounds on this processes are [|l2h 



We will carry out our calculations for the branching ratio in the mSUGRA scenario, where 
there are only four parameters that will fix the low energy values of the slepton masses 
Ml/2, iTT'O, ^0; tanP, sign{n) but our fit for the fermion masses was carried out for the tan/3 = 55 
so we will stay with this value. In the FIG|2lone can see dependence of the branching ratios on the 



5x 10" 



Ix 10" 



5x10" 



Ix 10" 



5x 10" 




Ml 



500 



1000 



1500 



2000 



FIG. 2: Branching ratios BR{n — > 67) for different values of mo (black 300 Gev, green 400 Gev, blue 500 
Gev), ^0 = 0, tan(3 = 55, sign(/i)=l 



Ml/2 for the fixed values of mo, tanjS, Aq, sign{fj.) . We note that branching ratio for /i — e + 7 for 
almost the entire parameter range of our model is above 10^^^ a value which is in the accessible 
range of the ongoing MEG experiment ! 1 2|] . 



VII. COMMENTS 



We add a few comments on the model described before closing: 

(i) In this model, the leading order proton decay operator is After GUT symmetry break- 

ing this leads to the effective strength A ~ -ttt-. Naively this is of order 2 x 10 , bigger than the 
present upper limit but is a considerable improvement in the naturalness. It could also be that the 
GUT vev could arise mainly from Ai with (^42) being an order of magnitude smaller. This would 
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then give the desired suppression to proton decay. In that case this will be the dominant graph for 
proton decay. Note that there are no Higgsino mediated diagrams for proton decay in this model. 
In addition, there is the gauge exchange diagram, present in all SO(IO) GUT models. 

(ii) The — e + 7 appears to be the only other low energy test of the model which is similar 
to such models. 

(iii) For the choice of parameters used in fermion mass fitting the neutrino mixing angles and 
mass differences could have any values in the allowed region. 

VIII. CONCLUSION 

In conclusion, we have presented the minimal SO(IO) 16-Higgs model for fermion masses 
where the problem of extreme fine tuning of higher dimensional Planck scale induced proton 
decay operators has been considerably ameliorated by the presence of discrete symmetries so that 
in the end, we only need to tune down the coupling only by a factor of 10~^. In this sense it is a 
more natural model We exhibited a fit to all fermion masses and mixings including neutrinos in 
this model to show that it can indeed be a realistic description of nature. 

This work is supported by the National Science Foundation Grant No. PHY-0652363. This 
work was presented at the GUT2007 workshop at Ritsumeikan University in Japan by R. N. M. 
and we like to thank the participants at this workshop and M. Ratz for comments. 
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